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Abstract: 

The electrophysiology of /9-cells is at the origin of insulin secretion. /9-cells exhibit a com- 
plex behaviour upon stimulation with glucose including repeated and uninterrupted bursting. 
Mathematical modelling is most suitable to improve knowledge about the function of various 
transmembrane currents provided the model is based on reliable data. This is the first attempt 
to build a mathematical model for the /5-cell-electrophysiology in a bottom-up approach which 
relies on single protein conductivity data. The results of previous whole-cell-based models are 
reconsidered. The full simulation including all prominent transmembrane proteins in /5-cells is 
used to provide a functional interpretation of their role in /3-cell-bursting and an updated vantage 
point of /3-cell-electrophysiology. As a result of a number of in silico knock-out- and block- 
experiments the novel model makes some unexpected predictions: Single-channel conductivity 
data imply that calcium-gated potassium currents are rather small. Thus, their role in burst in- 
terruption has to be revisited. An alternative role in high calcium level oscillations is proposed 
and an alternative burst interruption model is presented. It also turns out that sodium currents 
are more relevant than expected so far. Experiments are proposed to verify these predictions. 
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1 Introduction 



The electrophysiology of (3-celh is most relevant for understanding possible regulatory targets 
of insulin secretion. Exocytosis of insulin loaded granules is, among others, governed by in- 
tracellular calcium signalling. Calcium dynamics exhibits very specific patterns in reaction to 
increased glucose concentrations which includes repeated bursting of the membrane potential. 
Bursting denotes oscillations of action potential-like depolarisations, which are regularly inter- 
rupted and set in again after a phase of silence. Thus, oscillations occur on two distinct time 
scales: Fast oscillation of the membrane potential on the scale of 100 milliseconds and slow 
oscillations of intracellular calcium on the time scale of 10 seconds, i.e. in the rhythm of the 
repetition of bursting events. 

In a minimal model of oscillations of the membrane potential (1) it was found that an activation- 
delay between voltage-gated potassium and calcium channels gives rise to this bursting be- 
haviour. Sodium, the major player in neurons, is generally believed to be of minor importance 
(2). Therefore, modelling work was concentrating on potassium and calcium currents (3-5). 

Initiation of bursting events was found to be related to ATP-driven potassium channels (6, 7) 
which induce a relevant outflow of potassium under resting conditions. When glucose - and 
via glucose metabolism - ATP is increased, the channel conductivity is inhibited (8, 9), which 
leads to the initial depolarisation. Then voltage-dependent calcium currents are activated and 
the delayed response of voltage-dependent potassium currents leads to the bursting event. 

Simulations of the /?-cell-electrophysiology have shown how, in principle, repeated bursting can 
happen (1-3,5). And the found mechanisms are in agreement with experimental results. How- 
ever, none of the presently existing models fulfill the following requirements: (i) All important 
ions and membrane proteins are explicitly modeled including their activation and inactivation 
dynamics; (ii) the /3-cell exhibits a steady state and is stimulated by changes of glucose con- 
centrations; (iii) the dynamics of the membrane protein activity are fully derived from protein 
experiment. Only if all these requirements are respected, the exact role of the different mem- 
brane proteins might be disentangled. 

A major problem in developing such a complete simulation tool is related to the large variety 
of measured whole-cell conductivities in /3-cells. In fact the density of the membrane proteins 
is itself a dynamical quantity. Their dynamics vary for different cell-types so that the data used 
for modelling are restricted to measurements in /9-cells. Not all membrane proteins have been 
subject to quantitative measurements in /5-cells and the /?-cells under consideration can be in 
various states of protein expression. Thus, modelling work was restricted to fit the activity- 
dynamics of the membrane proteins to the behaviour which was expected. 

In order to overcome this problem, the present simulation exclusively relies on data of single 
protein activity. Such an approach was used in calcium modelling of neurons (10) before and 
is now applied for the first time to /5-cell electrophysiology. Today there is a rather complete 
knowledge about single proteins conductivity and opening dynamics available which allows a 
bottom-up approach starting from the single protein level. 

To this end the single protein activity dynamics is separated from the protein density. It is 
assumed that the single protein dynamics is cell-independent. In other words a voltage-gated 



2 



potassium channel has the same dynamical properties in a neuron and in a /3-cell. This assump- 
tion might be set in question if cooperation of membrane proteins with other cell structures 
would change the activity dynamics. However, the universality of single membrane protein 
properties can be considered to be a good approximation because the measured characteristics 
of the cells show relatively small variations despite them stemming from different cell-types. 

The conductivity and opening properties of single proteins are exactly implemented into the 
model and then are multiplied by the surface density of the respective proteins in /3-cells. The 
latter are not precisely known and are the free parameters of the simulation. In comparison 
to the huge parameter space in simulations relying on whole-cell conductivity measurements, 
this is a rather small set of unknown parameters. This improves the predictive power of the 
mathematical model. The behaviour of the /9-cell in terms of dynamics of membrane potential 
and ion-concentrations is emerging from the single protein level, thus, coupling the molecular 
to the cellular level. 

The framework of the simulation is presented in the methods section. The single membrane 
protein properties are collected in the supplementary material. These define the current dynam- 
ics as used in the simulation. The full model as introduced in Section |2] contains the model of 
a most prominent original work (3). The bursting behaviour as found in Sherman et al's work 
is reproduced as a special case of the full model and is used as test for the simulation. Then the 
full simulation is employed to set up a steady state of the /9-cell, and to allow for stimulation 
of the /?-cell with glucose. The protein densities are determined to find a realistic behaviour 
of the cell - where the single membrane protein properties are not touched. On the basis of 
this simulation various in silico experiments are performed and compared to real experiments if 
available. 

The in silico experiments include knock-out experiments and partial or total block of membrane 
proteins. Calcium-driven potassium channels are modeled according to data that are used for 
/5-cell- simulations for the first time. The measured dynamics leads to the conclusion that this 
channel is of minor importance at normal glucose levels, while it is most important for the (3- 
cell-behaviour at supra-large glucose levels. In this regime calcium-driven potassium channels 
drive uninterrupted bursting (7). LVA-channels, which have not been considered in other mod- 
els, NCX, and PMC A are suggested to be important for the interruption of bursting events at 
normal glucose levels. Also the importance of a dynamic reversal potential, which was mostly 
neglected in other work, is pointed out. It is found that the role of sodium currents was un- 
derestimated so far, and their impact on /5-cell behaviour is explained. The paper is concluded 
with a novel vantage point on the interplay of the different ion-conducting membrane proteins 
during repeated and uninterrupted bursting. Novel experiments are suggested to verify the new 
vantage point. 
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2 Methods 



2.1 Dynamics of ion concentrations and membrane potential 

Rate equations are used to describe the dynamics of intracellular ion concentrations. External 
ion concentrations are assumed constant. The ions treated explicitly are sodium N, potassium 
K and calcium C. The dynamics depends on the conductivity of various membrane proteins 
which each correspond to one term in the following differential equations. 

= — (PNa,V-^Na,V + 2pNa,K-^Na,KttNa,K + PNCX-^NCXttNCX + <^Na) 

dK _ ^ f ^ r r 4 ^ ^ 

— — I PK,ATP-'K,ATP + PK,V-'K,V + PK,Ca-'K,Ca — g «Na,KPNa,K-'Na,K + 

dC i f ^ ^1 

PCa,L-'Ca,L + PCa.TJCa.T — TT^^Ca^NCXPNCX-'NCX 



dt ZcaF (1 + Xc) V 3 

+PPMCA-^PMCA + Jca) ■ (1) 

Lj. denote the electrical single membrane protein currents, where x specifies the type of mem- 
brane protein as defined in Figure [T] (see also Tab.[T]in the supplement). Negative currents are 
defined to bring positive ions (charges) into the cell. As the currents are not defined as ion cur- 
rents but as electrical currents they have to be weighted by factors denoting the valence of the 
ion zca = 2 and the stoichiometry of the protein a^. The single membrane protein currents and 
their dynamics are discussed in the supplementary material Section 15.11 and are schematically 
visualised in Figure [B 

Px denote the corresponding surface densities of the membrane proteins. = 3 is a surface to 
volume factor that relates the current through the membrane to a change of the ion concentration 
in the intracellular volume. The cell geometry is assumed to be approximately spherical which 
implies ^ = 3/r where r is the cell radius. 

Note that a calcium buffer is included in the equation for calcium. The calcium buffer is as- 
sumed to bind and dissolve calcium faster than the typical time scale of the ion dynamics. Thus, 
the kinetic equation for the buffer can be solved in the steady state approximation (rapid buffer 
approximation). This implies an additional factor (1 + Xc) in the equation for calcium, where 

with Co the total concentrations of the buffering calcium binding sites. The most prominent 
buffer is calmodulin having more than one calcium binding site. Kc is the dissociation constant 
of the calcium binding sites. C is the free calcium concentration. The total calcium concentra- 
tion is given by 

The fraction of free calcium is fca = and adopts values between 0.1 and 1%. 
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Figure 1: Schematic representation of ion-conducting transmembrane proteins: The 
plasma membrane is visualised as horizontal grey line on which a set of transmem- 
brane proteins are shown. The shaded boxes show the acronyms of the proteins and 
the quantities that control their gating/activity in the model introduced below. The ions 
attributed to the arrows denote the direction of ion flow (assuming normal electrochem- 
ical gradients). The calcium buffers symbolise all l<inds of intracellular calcium binding 
sites, including calmodulin and organelles acting as calcium stores. The acronyms of 
the membrane proteins are also defined in Tab.[l| in the supplement, where also the 
activity dynamics are explained in detail. 



<^K,Na,Ca dcnotc the leakage currents. They are derived from the steady state limit (see Eq.[24lin 
the supplement) to guarantee a resting state of the cell. In the resting state all other currents are 
balanced by the leakage currents, so that the net current vanishes and the ion concentrations are 
constant. 

The ions build up the membrane potential V. All non-electroneutral currents introduced in Eq.[T] 
show up in the equation for V: 

- (PNa,K-^Na,K + PK,ATP-^K,ATP + PK,V-^K,V + PK,Ca-?^K,Ca + PNa,V-?^Na,V 
+PNCX-^NCX + PPMCA-^PMCA + PCa,L-^Ca,L + PCa,T-^Ca,T 

+ Jk + ^Na + Jca) • (4) 

No additional leakage current is needed because in the steady state the right hand side vanishes 
by construction. Cm is the membrane capacitance that relates the electrical currents to changes 
of the membrane potential. 



dV 
'dt 
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2.2 Reversal potentials 



Reversal potentials depend on the chemical gradient and the membrane potential over the cel- 
lular membrane. They cannot be considered to be constant (as assumed in most other models). 
As the reversal potentials linearly enter Ohm's law (which is used to approximate ion current 
through open pores) this has to be accounted for. The Nemst-equation is used to calculate the 
correct reversal potential during dynamical changes of the electrochemical constellation: 



where x = K, Na, Ca, and R = 8. 315 J /{Kmol) the Rydberg (molar) gas constant. The 
temperature T is assumed to be the physiological temperature. Xcxt is the external concentration 
of the ion (in thermic bath approximation). They are calculated from the intracellular ion- 
concentrations at rest state to fit the known reversal potential Vk = —75mV, ^Na = +80mV, 
V^ca = +l28mV. 

The straight forward application of Nemst equation, which is based on thermodynamical no- 
tions and only involves the electro-chemical gradient of the ion under consideration, thus, ne- 
glecting other ions and not treating diffusion of ions correctly, is not always justified. In the 
case of calcium a more exact theoretical approach based on the Goldman-Hodgkin-Katz equa- 
tion leads to a non-linear I-V-relationship (see e.g. (11)). This equation is derived from the 
theory of electro-diffusion through a membrane and involves the influence of other ions on the 
reversal potential. This more precise I-V-relationship is approximated by an Ohm's law like 
current, which is only justified for the linear piece of the curve. A good approximation is to use 
a corrected reversal potential of Vca = +50mV. Up to this rather large membrane potential 
the I-V-relationship is, indeed, approximately linear. However, one should keep in mind that 
for depolarisation beyond 50mV (thus of AV^ = 120mV relative to the resting potential) the 
calcium currents are not correctly represented. 

2.3 The basic model assumptions 

The following list summarises the most important features, assumptions, and limitations of the 
/3-cell-model: 

Universality of single membrane protein measurements: 

It is assumed that the measured single membrane protein conductivities and opening dynam- 
ics are valid irrespective of the cell-type. This considerably increases the data basis of the 
simulation. In view of consistent measurements of the same protein in different cells this ap- 
proximation is reasonable. These parameters are only varied in order to investigate robustness 
issues. 

Stability of the j3-cell-resting state: 

The leakage currents are assumed to equilibrate all other currents if all cell properties are at 
their resting values. This implies that the cell exhibits a stable steady state and that the cell 
returns to this resting state after stimulation. 




(5) 



6 



Glucose- and not ATP -mediated activation: 

Glucose is assumed to directly impact on the conductivity of the K,ATP-channel. Thus, the 
/3-cell can be stimulated by changes of the glucose level. However, the dynamics of ATP (6, 
12) which mediates the effect of glucose on K,ATP-channels (13) is not considered and left 
for future improvements of the simulation. Also any feedback between ATP and calcium is 
neglected (14). 

Averaging over the whole-cell: 

Inhomogeneities of the cell are averaged out. The ion-concentrations are assumed to be average 
quantities over the whole-cell. 

Neglect of calcium-induced-calcium-release: 

Calcium-induced-calcium release has been suggested to contribute to the calcium oscillations in 
response to glucose-induced potassium currents (15). However, blocking of related mechanisms 
does not affect observed oscillations in pancreatic /3-cells (16, 17). Oscillation in the stores 
might be orchestrated with oscillations in the cytosol (17). The present simulation neglects the 
effects of calcium-induced-calcium-release and assumes that the endoplasmatic reticulum acts 
mainly as an additional buffer. Thus, the binding site concentration cq is assumed to be larger 
than expected by calmodulin alone. 

Membrane channels follow Ohm's law: 

The ion flow through an open pore is assumed to follow Ohm's law which implies a linear 
voltage-current-relation. This is valid for most ion channels. Deviations are found only for 
extreme membrane potentials which are not considered within this investigation. 

The (3 -cell is isolated: 

In vivo /3-cells are in the context of islets and connected with each others by gap-junctions. This 
connection might influence the susceptibility of stimulation. Indeed, /5-cell- stimulation is more 
difficult in isolated cells (see e.g. (7)). The present simulation neglects effects stemming from 
being embedded into islets. 

3 Results 

3.1 Evaluate the Sherman et al. 1988 model 

The model of Sherman et al. (3) can be considered as a classical reference model. It is charac- 
terised by the following properties: 

• Membrane potential follows a classical ordinary differential equation with capacity 5.3 IpF 
and three ionic currents (no leakage and no sodium currents). 

• A voltage-gated potassium channel (corresponding to K,V here) is described by a sig- 
moidal function for its asymptotic value and a differential equation for the opening prob- 
ability (same as here). Inactivation of K,V-channels is not considered. 

• A calcium-gated potassium channel (corresponding to K,Ca here), which is described by 
a Hill function with Hill coefficient nx.ca = 1- 



7 



• A calcium channel (corresponding to Ca,L here) including inactivation and fast activation. 

• Constant reversal potentials are assumed (no Nemst equation is used). 

• The steady state of the system is defined as zero current of the three currents contributing 
to the voltage equation (thus no leakage currents are considered). 

• An equation for calcium ions including two terms and a constant factor / = 0.001 com- 
pensating for the buffer (thus, no buffer dynamics). One term accounts for Ca,L-currents 
the second for calcium removal (corresponding to PMCA here). 

• PMCA-activity is not described by a Hill-function but increases linearly with the calcium 
concentration. Note that in (3) the calcium extrusion term is assumed to be electroneutral, 
which is not the case as, typically, calcium is removed by PMCA. This pump is not elec- 
troneutral. Thus the PMCA-current also enters the equation for the membrane potential 
Eq.lHin the present model. 

The full model as introduced in Section[2]is reduced to these properties and the results in (3) are 
reproduced (data not shown). In a second step some parameters are modified to experimentally 
known values: capacity, resting concentrations, reversal potential of calcium. All these changed 
values are collected in Tab.[2l-[5]in the supplement. Surprisingly, all modifications together only 
moderately alter the results (see Figure [2l). 

The channels in the initial state of the simulation are not equilibrated (thus not in a steady 
state). This induces a certain inflow of calcium through Ca,L which depolarises the cell. As 
the K,V current reacts with a delay the potassium outwards current comes later and allows 
for an above threshold depolarisation. This leads to the oscillation on a higher voltage level. 
Now, free calcium accumulates in the cell due to repeated inflow. Note that the largest part of 
the incoming calcium is assumed to bind to the buffer, such that calcium accumulates much 
slower than according to the calcium inflow. More calcium induces increasing gK,Ca (note 
that the K,Ca-currents used in (3) are not in agreement with single protein experiments) and 
PMCA-activity, which in turn increases the potassium and calcium outflow, respectively. The 
oscillations are suppressed for some threshold K,Ca and PMCA conductivity. From here the 
whole process restarts. The essential properties are the delay between Ca,L- and K,V-current, 
the calcium buffering which induces a slow calcium increase, and a membrane protein which 
reacts to the increased calcium level by dropping oscillations. 

The back-reaction of PMCA activity on the membrane potential was neglected in (3). If this is 
included into Sherman et al's model, bursting disappears. The back reaction of PMCA activity 
on the membrane potentials induces hyperpolarisation of the cell. This infers that K,V and K,Ca 
(in view of the reversal potential of Vk = —75mV) but more importantly Ca,L (because of 
very low opening probability) are almost inactive. In this situation PMCA reduces intracellular 
calcium below its resting state while further hyperpolarising the cell. This is a problem related 
to lack of steady state in Sherman et al's model. The cell does not go back to resting state after 
a bursting process. In conclusion, neglect of a back-reaction of the process removing calcium 
from the cell on the membrane potential turns out not to be a justified approximation. 
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Figure 2: The simulation is consistent with Sherman et al's model: A simulation of the 
model of Sherman et al. (3) with adapted parameters. The membrane potential (left 
panel, grey line) exhibits repeated bursting events. Intracellular calcium concentration 
is plotted as ratio to the resting concentration of O.lyuM (left panel, black full line). It 
rises during bursting events and decreases during the silent phase. The contributions 
of the transmembrane currents are shown in the right panel (the inset lists the currents 
in the sequence of vertical appearance). Oscillations occur between voltage-gated 
potassium (K,V right panel, light grey line) and calcium channels (Ca,L right panel, 
darl< grey line). K,Ca- and PIVICA-activity follow the calcium level (right panel, medium 
grey and blacl^ line, respectively. 



3.2 /?-cell-bursting with calcium, potassium, and sodium currents 

The second simulation includes all membrane proteins as introduced in Section [2] and, thus, 
extends the model (3). The characteristics of activation and inactivation are taken from single 
protein experiments. The densities are adapted to find reasonable electrical bursting. The re- 
sulting set of densities is given in in the supplement Tab.[5](full model). First, the steady state 
is established and then the cell is stimulated by 10m Af glucose. The result is shown in Figure [3] 
and is in good agreement with corresponding measurements (7, 18). The whole-cell peak ion- 
currents are 12, 50, and 70pA for sodium, calcium, potassium, respectively, which is also found 
in experiments. The whole-cell K,ATP-conductivity is InS at ImM glucose (derived from 
Figure |4] as 47rr^pK,ATP5'K,ATp(l — 5'k,atp)). which compares to 1 — 3nS in experiment (19). 
An analysis of the contributions of the different membrane proteins induces a picture which 
in parts coincides with the one proposed by Sherman et al. (3). However, there are also some 
important modifications which will be analysed and isolated in subsequent knock-out and block 
experiments. 
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Figure 3: Repeated bursting including all transmembrane proteins: The full model for 
the electrophysiology of the /3-cell is used to simulate repeated bursting after activation 
with glucose. Before stimulation the /3-cell exhibits a stable steady state. Activation with 
lOmM glucose at t = 3s (left panel, dashed line) first leads to instabilities. Then the 
membrane potential (left panel, grey line) shows regular repeated bursting activity. The 
burst to silent ratio and the bursting frequency are realistic. The intracellular calcium 
level (left panel, black full line) increases and decreases in cycles corresponding to the 
bursting events. The whole-cell contributions of the different currents are shown in the 
right panel (the proteins in the inset are given in the sequence of current contributions). 
Note that here and in the corresponding subsequent figures the electrical currents, not 
ion-currents, are shown. Therefore the NCX-current appears with negative sign. 



3.2.1 High resolution of bursting events 

A single burst event in Figure [3] is analysed in more detail. To this end a high resolution figure 
of a burst event is depicted in Figure [5l The analysis is started when the membrane potential has 
reached its minimum. It can be clearly seen that the burst event is initialised by a K,ATP-current. 
The Ca,L-current quickly follows. The K,V-currents are suppressed in this phase of the burst 
event. This is not related to the delayed activation (which causes 50ms at the most) but to the 
general single protein activation properties. These do not allow for a strong K,V-current near the 
reversal potential of potassium. K,V-currents become relevant only after strong depolarisations 
which are induced by strong Ca,L-currents. The delay of K,V-response with respect to the Ca,L- 
dynamics leads to stable oscillations. Calcium levels increase with each depolarisation spike. 
This induces small K,Ca-currents, involving a slightly decreasing frequency of spikes during 
each burst event which is in agreement with experimental data (7, 18). 

Interruption is mediated by a concerted action of NCX, PMCA, K,Ca and Ca,T in response to 
the increasing calcium concentration. But on top of changed conductivity the reversal potential 
plays an important role. For calcium it changes from 22mV at the beginning of a burst to 
OmV at the end of a burst. This reduces the calcium currents through voltage-gated channels 
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Figure 4: Opening/activity states of the membrane proteins: The States of the transmem- 
brane proteins during repeated bursting as given in Figure [3] are shown as probability 
of activity or opening. The left panel collects all calcium-relevant proteins and the 
K,ATP-channel. The right panel all sodium-relevant and other potassium-conducting 
proteins (note the different scale in the right panel and the shifted open probability for 
K,V-channels). The inset lists the probabilities in the sequence of vertical appearance 
and gives a hint how the probabilities are calculated (see supplement Section [5711 for 
details). The simulation protocol underlying this simulation is explained in the figure 
legend of Figure [3l Note that the PIVICA-activity is at its limit in an actively bursting 
p-ceW. 



by about 35% (at V = —AOmV) and is a relevant contribution to the interruption of bursts. 
Indeed, replacing the Nernst-equation by constant reversal potentials (the ones corresponding 
to the resting state) leads to uninterrupted bursting (data not shown). 



3.3 Knock-out of K,Ca-protein 

The calcium-gated potassium channel K,Ca is for the first time implemented according to sin- 
gle membrane protein characteristics (20), which differ strongly from the values used in other 
models. Currently, no model for the K,Ca-opening dynamics exists. The measured opening 
dynamics depends on the membrane potential and the calcium concentration in a complicated 
way. A detailed analysis of the measured data (see supplement Section [5X7]) has revealed that, 
in addition, the asymptotic opening probability for a given calcium concentration also depends 
on the membrane potential (see Eq.[2T[). 

The novel model for K,Ca-opening probability is used for the /5-cell-electrophysiology and it 
is predicted that under normal conditions the K,Ca-currents remain rather small because of 
a low opening probability. This finding suggests that K,Ca might even be unnecessary for 
repeated bursting. K,Ca is believed to abort burst-events when calcium is sufficiently increased. 
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Figure 5: High-resolution view on a burst event: The third bursting event in Figure [3] is 
shown in higher resolution. All shown curves have the same meaning as explained in 
the legend of Figure [3l 



In principle, this can also be achieved by PMCA which also directly depends on the calcium 
level. K,Ca is switched off (pK,Ca = 0) and other densities have to be slightly adapted to 
Pk.atp = 0.13//im^ and Ppmca = 1420//^im^ in order to keep regular bursting (see Figure©. 
Note that this simulation does not correspond to a block of K,Ca. The /9-cell is constructed 
without K,Ca, thus, this simulation corresponds to a knock-out experiment. 
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Figure 6: /9-cell-bursting without K,Ca-channels: Simulation without K,Ca-channels and 
with protein densities adapted to pk.atp = O-lS/z/m^ and ppmca = 1420//im2. Bursting 
frequency and the reached calcium level are both higher compared to Figure [3l Lines 
are explained in the legend of Figure [3l 
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In the normal (3 -cell a late increase in K,Ca- activity happens for high calcium levels (at the end 
of each burst, see Figure [3] and Figure Sj) leading to potassium efflux. This counterbalances the 
calcium influx and drops depolarisation. In the present simulation without K,Ca, PMCA takes 
over this function. However, PMCA-activity is already at its limits during burst (see Figure H]). 
This is also found in the K,Ca-knock-out experiment (not shown). 

The K,Ca-knock-out experiment suggests that for low calcium concentrations PMCA is the 
main interrupter of bursts. The interpretation that K,Ca is important for burst interruption is only 
supported for calcium concentration above 500/iAf. In this regime PMCA is not sufficiently 
flexible to further increase its activity and, thus, to interrupt bursting. K,Ca is much more 
flexible in this respect. It will be seen later that at high calcium concentration K,Ca acquires 
other functionality. 

3.4 Stimulation by raising external potassium 

The effect of increased external potassium levels on /3-cells is investigated. Starting from a 
/5-cell in its resting state at normal external potassium A'cxt = 5.7mM, external potassium is 
increased to different values using a steep sigmoidal function. The /3-cell which had adapted to 
the lower potassium level is stimulated because the reduced chemical gradient pushes potassium 
into the cell and depolarises it. 

For moderately increased external potassium the cell first exhibits a smooth depolarising cal- 
cium current. Larger i^cxt then leads to K,ATP-induced action potential-like events involving 
potassium and calcium currents, and finally into a burst-like activity at A'cxt = 8mA'/. How- 
ever, this remains a unique event (data not shown). The burst is interrupted by K,Ca-activity 
at a rather high and constant calcium level around IfiM. Despite a new depolarised equilib- 
rium state which normally would induce action potentials, repeated bursting is prohibited by 
the constant strong activity of PMCA and K,Ca. 

3.5 Stepwise increase of glucose levels 

The full model is used to investigate the /5-cell behaviour after stimulation with increased glu- 
cose levels. This inhibits K,ATP-potassium currents and depolarises the cell. Depending on the 
target glucose concentration between the resting value 70 = ImM and 30mM this becomes 
the germ of bursting events. The /9-cell reaction is shown in Figure |7] for some examples. Up to 
5mM glucose only a smooth increase of the membrane potential is found which then stabilises 
on a higher level. A single action potential appears at a threshold of 6mM which is in agree- 
ment with experiment (18). The spikes then turn into regular firing with increasing frequency 
between 7 and 9mM. The high-frequency spikes then are interrupted and replaced by repeated 
bursting around lOmM. The frequency of repeated bursting, the burst to silent ratio, the burst- 
ing amplitude, and the bursting baseline at 7 = lOmM are within the range of experimental 
results (7, 18). For even larger glucose levels uninterrupted bursting is found which is also in 
agreement with experiment (7). 

As in islets (18) calcium levels constantly increase with increasing glucose levels (see Figure 
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Figure 7: The /5-cell reaction to increased glucose levels: Glucose Stimulation is increased 
in several steps 7stimuiation = 5,6,8, 10, 15,30mM (panels in reading order). The mem- 
brane potential (grey line) and the calcium concentration (black line) are shown. Re- 
peated bursting is found for stimulation with about lOmM glucose. Less glucose only 
induces single spil<es, more leads to uninterrupted bursting. 



|7]). This is not the case if the simpler model for K,Ca-activation is used (data not shown), which 
underlines the necessity of using the improved novel model for K,Ca-dynamics. 



3.5.1 Steepness of glucose increases 

In experiments it is observed that the initial spike or burst is more intense than the subsequently 
following ones - if any. According to the simulations this is related to the steepness of the switch 
to the higher glucose level. With smooth sigmoidal functions the initial spike can be fully sup- 
pressed. Similarly, long initial bursts can turn to a normal duration. However, the steepness of 
the glucose increase changes only the initial action potentials or pre-bursting events. The long- 
term bursting behaviour does exclusively depend on the asymptotic glucose level irrespective 
of how quick it is reached. 



3.5.2 The role of K,Ca-currents at supra-large glucose 

There is an unexpected diversity of model behaviour at high glucose levels in dependence of the 
model assumptions. Uninterrupted bursting for supra-large glucose levels is only found with the 
novel model for K,Ca-activation (see supplement Section [5.1.7|) . The simpler model for K,Ca- 
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activation as used in Figure [2l leads to suppression of any spike for large glucose, which is in 
contradiction to experiment. This underlines the necessity to rely on the measured electrophys- 
iology of single K,Ca-channels. A rather weak role of K,Ca-channels during normal electrical 
activity in inferred. K,Ca gets relevantly active only at very high calcium concentrations which 
are associated with high glucose concentrations. 

A more detailed analysis of Figure |7] reveals that uninterrupted bursting at supra-large glucose- 
levels is based on oscillations in the K,Ca-current. While at 7 = 12mM the K,V-current is 
dominant, at 30mM it is the K,Ca-channel which provides the largest contribution to the potas- 
sium current (see Figure [8]). Instead of K,V and Ca,L it is now K,Ca and Ca,L which are 
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Figure 8: K,Ca is the main oscillator at large glucose levels: The current contributions 
of K,Ca and K,V to uninterrupted bursting for glucose levels of 12 (left) and SOmM 
(right) are shown. While at 12mM the main contributors to the oscillations are K,V 
(light grey full line) and Ca,L (grey full line), these are K,Ca (black full line) and Ca,L at 
SOmM glucose. Note also that the contribution of K,ATP-channels is strongly reduced 
at supra-high glucose levels. 

swinging. This result suggests an additional role of K,Ca-channel in the electrophysiology of 
the (3-ce\\: It defines the electrical activity at supra-large glucose and calcium concentrations. 

3.6 Inhibit SERCA-activity 

Thapsigargin is a SERCA-inhibitor. Its effect is to elevate the intracellular baseline calcium 
concentration. It was found that thapsigargin can induce continuous firing in /3-cells (7). As the 
present simulation does not include the ER explicitly SERCA-block is mimicked by induction 
of its effect on the calcium baseline. To this end the calcium leakage current in a cell with 
established resting state is modified. 

A reduction of the leakage current to 20%, indeed, leads to an increased baseline calcium level. 
If then glucose is turned to 7 = lOmM continuous bursting is found. When the leakage current 
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is turned back to its resting state value normal repeated bursting is recovered (see Figure [9l). 
Continuous bursting is also found if the leakage current is reduced after repeated bursting being 
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Figure 9: Reduction of leakage current to mimic SERCA-block: The calcium leakage 
current is reduced to 20% at 3s to mimic SERCA-blocking (right panel, grey dotted 
line). This increases the calcium baseline (left panel, black line). At t = 20s glucose is 
turned to 7 = lOmM (left panel, black dotted line), which induces continuous bursting 
(left panel, grey line). At t = 80s the leakage current is turned back to normal and 
repeated bursting (as in Figure [3]) is recovered. At t = 110s glucose is turned back to 
its resting value. The right panel shows all but K,V whole cell potassium currents. Note 
that because of its stoichiometry Na,K-potassium current has negative sign and double 
absolute value with respect to the shown curve (right panel, black full line). 

already established at high glucose levels. Thus, the sequence of stimulation events is not 
important. 

From Figure |9] it can be read off that the K,Ca-current is increased during continuous bursting 
as compared to repeated bursting. In contrast, the K,ATP-current is smaller during continuous 
bursting. This is in agreement with the results found in (7), Figure 4 therein. However, the 
interpretation is different: The authors in (7) suspected that the reduced consumption of ATP 
(during SERCA block) leads to less inhibition of K,ATP-channels. The simulation suggests 
that the increased calcium baseline increases the necessary PMCA pumping activity in order to 
establish a steady-state. In fact, the consumed amount of ATP is more related to the absolute 
calcium level in the cell than to the number of blocked ATPases. Thus, the present simulation 
suggests that altered calcium concentrations and membrane potentials during block of SERCA 
reduce K,ATP-currents independent of the intracellular ATP-concentration. 
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3.7 Partial blocking of PMCA and NCX 



A very similar result (to Figure [9l) is found when the PMCA-current is reduced to 90%. The 
calcium baseline is increased and at lOmM glucose uninterrupted bursting is found. Substan- 
tially stronger inhibition of PMCA- activity suppresses bursting after an initial bursting event 
(data not shown). 

A reduction of the NCX-current to 80% also leads to uninterrupted bursting. However, the 
characteristics are different. As before the calcium baseline is increased. But, instead of an 
increased calcium level during uninterrupted bursting, the calcium level is decreased (data not 
shown). This points to two different mechanisms both leading to uninterrupted bursting. 

3.8 Block K,ATP with tolbutamide 

Block of K,ATP-currents was investigated in patch-clamp experiments with /5-cells in intact 
islets (7). In silico the application of tolbutamide corresponds to a reduction of the K,ATP 
protein density pk,atp in a cell that has established its steady state in the presence of normal 
K,ATP-presence. A reduction of K,ATP-currents to 75% leads to a result almost identical to 
Figure |9l An increased calcium baseline is established. This happens because the reduced 
potassium outflow depolarises the cell which, in turn, partially opens voltage-dependent calcium 
channels. The new equilibrium state is, thus, established on a higher level of calcium influx. 

When glucose is increased to lOmM uninterrupted bursting is induced. As in experiment (7) 
this is reversible: After release of partial K,ATP-block normal repeated bursting behaviour is 
reestablished. This happens in the simulation without any additional stimulation. In (7) the 
/3-cell was depolarised several times in voltage clamp mode. 

Note that an overexpression of K,ATP-channels also leads to uninterrupted bursting despite the 
fact that the calcium baseline is decreased (data not shown). Again this points to at least two 
mechanisms for uninterrupted bursting. For even more increased K,ATP-currents regular single 
firing is observed. 

Interestingly, a total block of K,ATP-channels immediately leads to continuous bursting even 
at the resting glucose concentration 70 = IniM (see Figure [TOl) . This supports that these 
channels trigger bursting events by inhibition of potassium outwards currents. The observed 
independence of the increased glucose level at t = 20s is a direct consequence of the fact that 
in the model K,ATP-channels are the only point of effect of glucose. ATP-ADP ratios are not 
explicitly monitored which would also change ATPase activities. 

It is possible to stimulate /?-cells by partial block of K,ATP-currents. Repeated bursting is 
induced at resting glucose level (see Figure [lOl right panel, for a 50%-block). When glucose is 
increased at t = 20s repeated bursting turns into uninterrupted bursting. 

3.9 Modulation of voltage-gated calcium channels 

While overexpression of LVA calcium channels Ca,T leads to similar results as shown in Fig- 
ure m block or knock-out of LVA channels has only minor effect on the bursting behaviour. 
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Figure 10: Total and partial block of K,ATP: K,ATP is fully blocked at t = 3s (left panel). 
This induces strong depolarisation, calcium inflow and continuous bursting. Att = 20s 
glucose is turned to 7 = lOmM, which does not influence bursting for lack of K,ATP. 
Att = 80s the block of K,ATP is released and repeated bursting is recovered (as found 
in Figure [3). Att = 110s glucose is turned back to its resting value. The right panel 
shows the result of partial block of K,ATP-channels to 50% using the otherwise same 
protocol. High frequency repeated bursting is found at resting glucose levels. 



Inhibition only weakly modifies the calcium baseline before stimulation with glucose. The am- 
plitude of spikes of the membrane potential is slightly increased and the reached calcium level 
is slightly shifted (data not shown). Thus, according to the present simulations, LVA have to be 
considered unimportant for /3-cell-bursting. 

However, this is a weak statement because for a different set of densities (with a larger fraction 
of LVA channels) the role of LVA calcium channels is more pronounced. In particular, the 
inactivation of Ca,T (note that HVA Ca,L-type calcium channels inactivate only slowly) might 
play a role for interruption of bursting events. In these simulations uninterrupted bursting is 
found for weak and for total inhibition of LVA-currents. The present simulation does not yet 
allow to make precise statements on the importance of LVA channels and on their inactivation. 
However, in can be stated that if Ca,T-type channels were expressed to a relevant amount in 
/5-cell, their inactivation properties would have impact on the interruption of bursting events. 
This is further supported by the observation that 70% of the calcium current in whole-/3-cell- 
experiments is inactivated after 100ms (19), which is impossible on the basis of Ca,L-channels 
that show only little or no inactivation. It would be most interesting to get data about the ratio 
of LVA and HVA calcium channels in /?-cells. 

Moderate inhibited or strengthened Ca,L-currents also lead to uninterrupted bursting combined 
with an increased or lowered calcium baseline, respectively (data not shown). For stronger 
inhibition (20% remaining) spiking activity disappears, which is consistent with the dominant 
role of the Ca,L-current for /3-cell-oscillations. 
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3.10 Modulation of sodium currents 



The sodium-potassium exchanger Na,K uses ATP to transport both sodium and potassium 
against their electrochemical gradient. Even though this membrane protein is only weakly ac- 
tive in the resting state, shows little dynamics upon /3-cell stimulation (variation of intracellular 
potassium concentration is relatively small), and induces only small currents, the bursting be- 
haviour is sensitive to partial blocking of the exchanger. Already at 70% of the Na,K-current 
uninterrupted bursting is found. This comes with a comparably small increase in the calcium 
baseline. Uninterrupted bursting is also found for increased Na,K-currents (130% or more). In- 
terestingly, the calcium baseline is decreased in this scenario, and, thus, cannot be at the origin 
of uninterrupted bursting. Note that due to the weak dynamics of Na,K-currents corresponding 
effects are not found in Na,K-knock-out experiments. 

Similarly the voltage-dependent sodium channel Na,V is generally not considered to play a ma- 
jor role in /3-cell electrophysiology. At a partial block to 80% (or lower) uninterrupted bursting 
is found together with a slightly reduced calcium baseline (see Figure [TT] left panel). Uninter- 
rupted bursting is also found if the Na,V-currents are increased by a factor 1.3 or higher together 
with a slightly increased calcium baseline (see Figure [TT] right panel). 
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Figure 11: Modulating Na,V-currents: Na,V-currents are inhibited to 90% (left panel) 
and increased to 130% (right panel) using the same protocol as in Figures [9] and [TOl 
The calcium baseline (full blacl< line) is reduced and increased, respectively. At lOmM 
glucose (see t > 20s) uninterrupted bursting is found in both cases with differing char- 
acteristics. While uninterrupted bursting with increased calcium baseline is similar to 
Figure H amplitude and burst baseline are different when coming together with a re- 
duced calcium baseline. This points to different mechanisms leading to uninterrupted 
bursting. 



Why is uninterrupted bursting found with reduced calcium baseline when either Na,K-currents 
are increased or Na,V-currents are reduced? Both, increased Na,K- and reduced Na,V-currents 
polarise the cell because intracellular sodium levels are reduced. This reduces the calcium 
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influx through vohage-gated calcium channels, and, in turn, reduces the calcium baseline. In 
the same way intracellular potassium levels are increased because of less potassium efflux. This 
corresponds to a pre- stimulation in complete analogy to the inhibition by K,ATP-channels. 

Thus, we have to distinguish low and high calcium uninterrupted bursting. /5-cells become 
prone to oscillations either by increased calcium baseline or by increased potassium baseline. 
The equilibrium of both is regulated by the sodium-currents. Despite the sodium conducting 
proteins playing an unspectacular role in /5-cell activity they nevertheless influence the electro- 
physiology in an unexpected strength. 

3.11 Modulation of potassium currents 

Inhibition of K,Ca-currents reduces the burst to silence ratio and the bursting frequency. Even 
at currents of 10% repeated bursting is still observed. Only when K,Ca-currents are totally 
blocked, bursting is fully suppressed. Increased K,Ca-currents lead to increased burst to silence 
ratios and at 150% to uninterrupted bursting (the characteristics correspond to those as found in 
Figure [m left panel). The calcium baseline and the membrane potential are not visibly altered 
by modulation of the K,Ca-currents (data not shown). 

Increased K,V-currents lead to uninterrupted bursting. A reduction of K,V-currents reduces the 
burst to silence ratio. Surprisingly, even with totally blocked K,V-currents, oscillations are still 
found with almost unaltered calcium baseline. The potassium current in the action potentials 
is then taken over by the K,Ca-currents. This shows that K,Ca- proteins have the potential to 
drive the oscillations during burst events - similar to the scenario found for supra-large glucose 
levels. 

3.12 Long-term stimulation 

Stimulation by increased glucose levels not only induce repeated bursting, the /9-cell is also 
shifted out of its resting state. If stimulation by glucose is applied for long durations, and 
assuming in a thought experiment all other cell properties unaltered during that period, the /3- 
cell might find a novel stable resting state which is in equilibrium with the high glucose level. 
The simulation in Figure [3]is repeated with lOmM glucose applied for one hour. Due to lim- 
ited responsiveness of Na,K-exchangers each burst event (involving Na,V-current oscillations) 
increases the intracellular sodium level. After half an hour sodium reaches a novel resting 
state at 170% of its resting concentration adopted to the ongoing stimulation (data not shown). 
Similarly the intracellular potassium concentration decreases to a new equilibrium state. Now, 
Na,K-activity and sodium inflow at each bursting event get into equilibrium. Repeated bursting 
is still present at this higher sodium-state. Only the baseline of bursts is increased by about 
20mV. 
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4 Discussion 



A bottom-up approach for the (3 -cell with quantitative ambitions 

A simulation tool for the electrophysiology of /5-cell was introduced. The model is strictly 
based on experimental data for single membrane protein activity or conductivity measurements. 
This is considered to be a first step towards a quantitative modelling of /3-cell electrical activity. 
Previous modelling work either structurally considered only those elements necessary to induce 
repeated bursting or omitted some essential protein properties. All previous models relied on 
whole-cell experiments. Thus, the present model can be considered to be the first bottom-up 
approach of the /5-cell starting at the molecular level and predicting emerging properties on the 
cellular scale. 

Improvements with respect to other modelling approaches 

The present model, besides on relying on single protein characteristics, includes several /5-cell- 
properties that have been neglected in previous modelling approaches. It was found that the 
back-reaction of the PMCA-currents on the membrane potential are very important and strongly 
alter bursting behaviour as found in (1,3). The inactivation properties (neglected in (2)) of the 
channels turn out to be relevant for burst interruption and the stability of oscillation. Finally, the 
Nemst-equation was used to dynamically calculate the reversal potential (which was assumed 
constant in (1, 3, 5)) and turned out to substantially change the currents during burst events. 

The endoplasmatic reticulum might change the dynamics 

Even though a block of processes related to the endoplasmatic reticulum does not affect /5- 
cell oscillations (17), calcium-induced-calcium-release might well modify the dynamics of rais- 
ing calcium. These modifications might lead to adaptations of protein densities. Also the sim- 
ple minded cell geometry has impact on the protein densities. Thus, even though the resulting 
whole-cell currents are realistic, the protein densities used in the simulations (see Tab.[5]in the 
supplement) are not considered as a quantitative prediction in the present stage of the model. 
Note that calcium pumps in the endoplasmatic reticulum will in first approximation act as an 
additional sink of free calcium which, in the present model, is effectively included in a large 
concentration of calcium binding sites. 

Characteristics of glucose-dependent (3-cell-stimulation are reproduced 
Despite the necessity for further improvement, the agreement with current experimental data is 
rather convincing and strengthens the relevance of the found results. The calcium response is of 
similar strength, and whole cell currents are in the right range (19). Bursting is found in the same 
regime of glucose levels as in experiment (7, 18). The glucose threshold for /3-cell electrical 
response is quantitatively reproduced to be around 5mM (18). The bursting frequency is in 
the range of experimental data. Uninterrupted bursting is found for large glucose levels as in 
experiment (7). Thus, it can be stated that the simulation respects the relevant characteristics of 
repeated bursting in /9-cells. 

The range of glucose leading to bursting is too small 

A more detailed analysis reveals quantitative deviations between simulation and experiment. 
The most relevant one concerns the range of glucose for which repeated bursting is found. 
While the simulation predicts 9 to llmM in experiment the lower limit is at 7mM (18) and 
the upper limit at 15mM (7). This discrepancy might be explained by a structural difference 
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between the experimental set-up and the simulation: The simulation models single /9-cells, 
while both cited experiments investigate /5-cell in a context of intact islets. Repeated bursting is 
difficult to be induced in isolated /3-cells, which can explain why the range of glucose is so small 
in the simulations. It has to be evaluated in future research whether the gap-junction-mediated 
contact to other /5-cells can enlarge the glucose regime for repeated bursting. 

A novel model of 13 -cell-bursting 

The fact that the main properties of /9-cell-electrophysiology are correctly described by the 
simulation encourages to draw a novel picture of how /9-cell-bursting emerges. This picture is 
presented as a sequence of events: 

• Glucose inhibits K,ATP: Increased glucose only gives a short pulse of inhibition of 
K,ATP-currents. This slightly depolarises the /9-cell which is further amplified by induced 
sodium and calcium entry through voltage-gated channels (Na,V; Ca,L; and Ca,T). 

• KjATP-currents increase again: Even though the opening of K,ATP-channels is inhib- 
ited the current increases upon cell depolarisation. According to Eq.[l4](see supplement) 
inhibition is over-compensated by the increased electrochemical gradient. 

• Ca,L and Na,V cause electrical spiking: If the stimulus is over-critical a spike of cell 
depolarisation is induced by voltage-dependent calcium and sodium currents. 

• Oscillations are induced by delayed K,V-response: The potassium current repolarises 
the cell, but this happens only with a delay tk^v It is well known that delayed response 
can lead to stable fast oscillations. 

• Calcium rises during the burst: During a bursting event calcium and sodium is con- 
stantly entering the cell while potassium is leaving it. While the effect on potassium and 
sodium is small, calcium rises because of a strong electrochemical gradient and the small 
resting concentration. The speed of calcium increase is controlled by the intracellular 
calcium buffer, and is much slower than the oscillations. 

• The burst is interrupted by a concerted action of Nernst-reversal potential, PMCA, 
NCX, K,Ca and Ca,T: During bursts calcium levels increase and the calcium reversal 
potential is reduced. The latter inhibits the calcium-currents, while the first increase re- 
polarising PMCA-, NCX-, and K,Ca-currents, where K,Ca-currents play a minor role at 
normal calcium levels (20). At a critical calcium concentration and reversal potential 
oscillations are suppressed. The interruption of bursts is facilitated by the inactivation 
dynamics of LVA calcium channels. 

• Uninterrupted bursting is driven by K,Ca: As also K,Ca acts with a certain delay this 
membrane protein has totally analogous potential to drive oscillations with Na,V- and 
Ca,L-mediated depolarising currents. According to the measured activity characteristics 
(20) oscillations are possible at large calcium concentrations only and can even drive 
uninterrupted bursting. 

• Repolarisation after burst events is governed by PMCA, NCX, and Na,K: In com- 
plete analogy to neuronal repolarisation, ATPases and exchangers are needed to reestab- 
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lish the resting concentrations. However, sodium and potassium are continuously chang- 
ing their level during long stimulation. 

Note that this model differs from older models in several aspects which are discussed in more 
detail in the following. 

Sodium currents are important regulators of [3-cell-electrophysiology 

Na,V-currents directly participate in spikes and bursting events. The activity of the /3-cell is 
responsive to all properties of the Na,V-channel. Thereby the half activation potential \4v[a,v 
has the greatest influence. Na,V-currents facilitate onset of action potentials due to their fast 
dynamics (the reaction time is the fastest of all considered channels), participate in the interplay 
with delayed potassium currents to induce oscillations, and also facilitates the end of action 
potentials by relatively fast inactivation. The interpretation of a tightly controlled interplay of 
Na,V with other membrane proteins is supported by the observation that uninterrupted bursting 
is observed for partially blocked and overexpressed Na,V-proteins. The latter bursting is char- 
acterised by high calcium, the former by high potassium. It turns out that Na,V-currents regulate 
the equilibrium between both ions in /5-cell-stimulation. According to the model results sodium 
channels are more relevant for the /5-cell electrophysiology than previously expected. 

Long-term stimulation gives rise to a third time scale 

The in silico experiment of long term stimulation of /9-cells with high levels of glucose has 
revealed three time scales of /9-cell-electrophysiology: 

1 . Delay of K,V or K,Ca currents with respect to calcium and sodium induced depolarisa- 
tions (10 - 100ms), 

2. slow increase in intracellular calcium level due to intense buffering which keeps the frac- 
tion of free calcium low (10s), and 

3. even slower modifications of intracellular sodium and potassium levels (10mm). 

The third time scale should be observable in experiment. 
Glucose only induces a tiny inhibition of K, ATP -currents 

In the model the only effect of increased glucose levels is inhibition of K,ATP-currents. How- 
ever, the inhibition turns out to be a tiny effect which quickly is hidden by the increasd K,ATP- 
currents in response to cell depolarisation. This opens the possibility of alternative effects of 
glucose via ATP-metabolism. Several channels and transporters rely on ATP-concentrations and 
may alter their activity in response to changed glucose levels. Some might well initiate bursting 
as already established for K,ATP. This might explain the recently observed insulin secretion and 
calcium dynamis in islets of mice lacking K,ATP-activity (21). 

Bursts are interrupted by NCX, PMCA, Ca, T and the dynamic reversal potential 
Interruption of bursting only weakly rely on K,Ca-currents because of small activity at phys- 
iological calcium levels. Instead, interruption is strongly related to NCX and PMCA activity. 
Mostly neglected in former research is the impact of a dynamic calcium reversal potential (see 
Eq.[5]) which turns out to change the calcium currents by about 35% during a burst. The inac- 
tivation dynamics of LVA calcium channels Ca,T also supports burst interruption. The latter 
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statement depends on the ratio of LVA- and HVA-channels in /9-cell, which would be important 
to be measured. In summary, burst interruption seems not to rely on a single type of membrane 
proteins but on a concerted action of various factors. 

Calcium driven potassium channels K, Ca can drive oscillations 

For the first time it is postulated that K,Ca is an active part of the bursting behaviour as well. 
Thereby the unknown K,Ca-time scale of activation is relevant: If the K,Ca adapts faster to 
the membrane potential (thus using a smaller rK,ca) the currents get out of phase and uninter- 
rupted bursting is destroyed. Therefore, the measurement of rK,ca is most important to learn 
about this predicted functionality of K,Ca-currents. According to the simulation results K,Ca 
acquire a major role for very large calcium concentrations where it takes over the role of K,V 
in oscillations. 

Falsification of the predictions of the simulation 

On one hand an improved version of the simulation model has to be developed in order to 
see whether the predictions of the simulations and the attributed roles of various membrane 
proteins for /3-cell electrical activity are robust. The improved model shall include the impact of 
intracellular calcium stores and explicitly model the dynamics of ATP. On the other hand some 
of the in silico experiments were not executed in real systems. The simulation tool might serve 
as a guide for the design of very focused and conclusive experiments. This is not restricted to 
the proposed experiments. On the contrary the author welcomes the challenge of the simulation 
model with further suggestions for experimental set-ups. 
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5 Supplementary material 



5.1 Conductivity of membrane proteins 



The (electrical charge) currents entering Eq.[T] and Eq.|4] are themselves dynamic quantities. 
Channels are gated in dependence on the membrane potential, on the presence of ATP, and on 
glucose- and ion- concentration. The resulting current depends on the electrochemical gradient 
over the membrane. So every membrane protein has its own characteristics and the conductivity 
of passive or active ion flow depends on different other parameters. In the following the different 
membrane proteins (as listed in Tab.[T]) and their gating/activity properties are described in more 



Na,K 

Na,V 

NCX 

PMCA 

K,ATP 

K,V 

K,Ca 

Ca,L 

Ca,T 



sodium-potassium exchanger (Na+ outwards, K"*" inwards) 
voltage-gated sodium channels 

sodium-calcium exchanger (Na+ inwards, Ca^+ outwards) 
Plasma-membrane-calcium- ATPase (Ca^"*" outwards) 
ATP-driven (and glucose-dependent) potassium channels 
voltage-gated potassium channels (also called delayed rectifier) 
voltage- and calcium-gated potassium channel 
L-type voltage-dependent calcium channels 
T-type voltage-dependent calcium channels 



Table 1: Acronyms of transmembrane proteins: Acronyms for the different plasma mem- 
brane proteins contributing to ion flow through the membrane are defined. Inwards and 
outwards specifies in which direction the ion flows, i.e. in or out of the cell. 

detail. Currents through pores are generally approximated by a linear voltage-current relation, 
which infers currents proportional to the electrical potential gradient (terms with V — V with 
V the reversal potential of the ion under consideration). For some channels and for extreme 
depolarisations there might be deviations from this approximation. 

The activation dynamics of active carrier proteins relies on the Hill-function 



H{x,Xh,n) 



X 



+ xl 



(6) 



which can be derived from chemical kinetics, x denotes some concentration, Xh is the corre- 
sponding concentration of half activation, and n is the Hill-coefficient. Large n make the switch 
between inactive and active states steeper. 

The asymptotic state of activation and inactivation is described by two functions of sigmoidal 
shape 



C^inact {x,Xh, A) 



1 



1 + exp{(xh - x)/k} 
1 

1 + exp{(x — Xh)/A} 



(7) 



Such sigmoidal functions are frequently used in literature and the parameters of these functions 
are determined in experiments. Xh is the concentration of half-maximum value and k and A 
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regulate the activation and inactivation steepness, respectively. Inactivation is a property of ion 
channels which will be handled in complete analogy to activation. In a first approximation 
ATPases driving ions against their chemical gradient are inactivated. Only the channels Na,V; 
K,V; K,Ca; Ca,L; and Ca,T (see Tab.[3]and|4l) are subject to inactivation. 

5.1.1 The sodium-potassium exchanger Na,K 

The sodium potassium exchanger transports sodium out of the cell and potassium into the cell. 
Both transport processes are against the chemical gradient, so that this process has to consume 
energy which is provided by ATP. The current is modelled as 

^Na.K = ^Na.K (1 " H{K, i^Na.K, J^Na.x)) H{N, i^Na.K, ^Na.x) , (8) 

where /Na,K denotes the maximum current. The quantities used for the Hill-function are ex- 
plained after Eq.[6l The first term describes a decreasing activity with increasing potassium 
concentration. The second one an increasing activity with increasing sodium concentration 
(see (22)). This approach follows approximations frequently adopted in other models: The 
adaption of Na,K-activity to sodium and potassium levels is assumed quick such that the steady 
state approximation is justified. A direct impact of glucose on Na,K-activity in /3-cells (23) and 
a voltage-dependence of Na,K-activity (22) is neglected. As ATP-dynamics are not modelled, 
the ATP-level does not influence Na,K-activity. 

Activation dynamics: The used parameters follow those found in (22) to fit experimental ac- 
tivation dynamics. Good agreement is found with riNa.K = 2 and -ft'Na.K = 33mA'/ as well as 
nNa.,K = 2 and -ft'Na.K = 20mM for dependence on potassium and sodium, respectively. 

Stoichiometry: The stoichiometry of the exchanger is 3:2, thus 3 sodium ions are exchanged 
against 2 potassium ions entering the cell, i.e. ctNa.K = 1-5 in Eqs.[T]and|4l As both ions have the 
valence 1 this corresponds to 1 positive charges leaving the cell for 2 potassium ions entering 
the cell. Thus /Na,K > enters the different equation with different factors: +1 for Eq.lH —2 
for K, and +3 for N Eq.E 

Turnover rate: The turn over rate of one ATPase measured as the consumption of ATP- 
molecules per second is /Na,K = 200ATP/sec ^ 0.00003pA in canine cardiac neurons (24), 
where the last equality assumes one exchange of 3 sodium and 2 potassium per ATP. 

5.1.2 The voltage gated sodium channel Na,V 

The voltage gated sodium channel opens when the cell is depolarised. The sodium ions then 
follow their electrochemical gradient to flow into the cell. The current is therefore composed 
of two factors, one describing the open-probability of the channel, the other describing the flow 
through the open channel: 

-^Na,V = ^Na,V5'Na,V5'Na,V {V — ^Na) ■ (9) 

^Na is the sodium reversal potential, and gNa.,v is the maximum conductivity. The open and 
inactivation probabilities g, /iNa,v are assumed to exponentially reach their asymptotic values 
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and, thus, follow the kinetic equations 




dt 



<^act{V, Vka.V, ^Na,v) ~ ^'Na.V 
TNa,V 

C"inact 

(V, H^Na.V, ANa,v) — ^Na,V 
6'Na,V 



(10) 



The parameters are explained after Eq.|7l It is TNa.v that defines the time constant of adaptation 
of the sodium channel (the speed of conformity change) in response to a changed membrane 
potential, while 6'Na,v determines the time scale of channel inactivation. 

Conductivity: The single channel conductivity of the voltage-gated sodium channel is (^Na.v = 
lApS as found in the squid axon (25). This value is consistent with findings in other systems. 
Taking the results in mice /?-cells (26), one has a single channel current of IpA at —30mV 
which (using VNa = 80mV) corresponds to 9pS. 

Activation properties: Activation is observed between —QOmV and —lOmV (26). Thus, a 
realistic approximation for the activation dynamics is VNa,v ~ —35mV and AtNa.v ~ 8mV. 

Inactivation properties: The role of voltage-gated sodium channels in insulin secretion and 
in membrane depolarisation is controversial. Depending on species, the half-maximal inacti- 
vation is found at —109mV in mouse /?-cells (26) (where 2.6mM calcium were used, and at 
less calcium the inactivation point was shifted to even lower voltage) and total inactivation is 
found at —AOmV (27,28). Thus, one may estimate PVNa,v ~ — lOOmV" with a steepness of 
ANa,v ~ 20mV (in order to get the observed total suppression of inactivation at —IbOmV, and 
total inactivation at —AOmV). This very low half inactivation potential if fully consistent with 
measurement in islet-/9-cells (19) and suggests an unprobable major role in the electrophysiol- 
ogy of /9-cells (26), which is not confirmed in the present simulations. 

Time scales: According to the sodium currents as observed in squid giant axons and to the clas- 
sical description by Hodgkin-Huxley the time constant of activation depends on the membrane 
potential (11). This is represented by 



with c = 11.5ms, a = AOmV, b = 50mV, and Vq = —70mV. Similar rNa,v are found in 
islet-/9-cells (19). Note that the same simulation results are found with the constant average 
value of TNa.v = 3ms. Inactivation happens on a time scale of 6'Na,v = 4.6ms which was found 
by two pulse experiments (11). 

Calcium permeability of sodium channels: Sodium channels in neurons are in general perme- 
able for calcium ions. Therefore, Na,V-channels are explicitly blocked with TTX in experiments 
to investigate the opening dynamics of voltage-gated calcium channels in order to separate these 
from the currents through sodium channels (see e.g. (29)). Calcium permeability of Na,V is not 
considered in the present approach. 



c 



(11) 



TNa,V — 



exp{(y - Vo)/a} + exp{{Vo - V)/b} 
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5.1.3 The sodium-calcium exchanger NCX 



The sodium-calcium exchanger takes advantage of the electrochemical gradient of sodium 
(pushing sodium into the cell) to transport calcium ions against its electrochemical gradient. 
This exchanger is important to reestablish the resting state calcium concentration after an ex- 
citable calcium influx into the cell. As in the case of the sodium-potassium exchanger a fast 
adaption of the activity to the calcium concentration is assumed. The current is derived from a 
corresponding Hill-function: 

-^Ncx = IncxH {C , C-!<(cx, n^cx) ■ (12) 

Typical values for the single membrane protein parameters can be found in (10): /ncx = 
-0m05pA (30) and Cncx = 0.0022mM (31). The HiU coefficient is assumed to be n^cx = 1- 
This is justified by the linear relationship between the observed current and the calcium con- 
centration as observed for low calcium concentrations in cardiac cells (32). However, the NCX- 
Hill-coefficient is controversial (2,33). 

Stoichiometry: The stoichiometry of NCX is 3:1, i.e. oncx = 3 in most tissues (31), thus 
three sodium ions have to enter the cell in order to expel one calcium ion. The net current of 
NCX is one positive charge getting into the cell for every expelled calcium ion. This infers that 
-^Ncx < 0. /ncx enters the different equations with factor +1 for the equation for the membrane 
potential Eq.H a factor —1 in the calcium equation, and a factor ttNcx in the sodium equation 
Eq.m 

Impact of sodium-concentration: Note that in the present description Eq.[T2]the sodium con- 
centration does not enter, even though NCX activity also depends on the sodium concentration 
gradient. Even though the sodium dynamics is modelled, it is assumed that the sodium concen- 
tration gradient is approximately constant. This turns out to be a reasonable approximation. In 
contrast, to neurons, in /9-cells peaks in the potential are dominantly generated by calcium and 
not by sodium. 

Impact of ATP: A dependence of the NCX activity on intracellular ATP concentration has 
been observed (34). This might be included in improved versions of this model that treat ATP 
dynamically. 

5.1.4 PMCA 

PMCA is an ATP-driven calcium pump extruding calcium from the cell to the extracellular 
medium, which has been characterised in /9-cells (35). In a first attempt the dependence on the 
ATP concentration is ignored and the ATP-concentration is assumed to be large enough to make 
the pump work optimally. Then the activity is mainly dependent on the calcium concentration 
in the cell. A suitable modelling approach is 

/pMCA = /pMCA-f/lC*, C'PMCA, 'T-PMCa) • (13) 

The maximum current /pmca is positive, as it carries calcium out of the cell. A delayed activa- 
tion of PMCA by increased calcium has been observed (36) which is neglected in the present 
model (using the steady state approximation). 
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Stoichiometry: One calcium ion is transported per used ATP molecule leading to the stoi- 
chiometry of opmca = 1 (37). This is surprising in view of the Hill-coefficient which was 
determined to be npMCA = 2 (38), which would suggest a corresponding stoichiometry (see 
also the discussion in (39)). 

Turnover rate: The turnover rate of single PMCA is in the order of 30 Hertz (30) corresponding 
to an activity rate of ka = 0.03/ms (this value was also used in (3)). This turnover rate can be 
translated into an electrical current: Every pumping event corresponds to a flow of 2 electrical 
charges e leading to Jpmca = = 60 ■ 1.6 ■ 10~^^C/s ^ IQ-^^A = lO'^pA. 

Half activation concentration: Typical values for the half activation calcium concentration 
are Cpmca = O.l/iM (see e.g. (40) Figure 3). Note that the values of half activation are 
found to depend on the calmodulin concentration and the isoform under consideration (40) 
which is neglected in the present approach. The 4b-isoform has a slightly larger half activation 
concentration of 0.16/uM (40). 

ATP consumption: If modelling ATP concentration dynamics PMCA has to be considered as 
a major consumer of ATP. Note that up to 40% of the ATP consumed to pump calcium out of 
the cell is used by SERCA-ATPase, thus, not in the plasma membrane but at the endoplasmatic 
reticulum (14). 



5.1.5 ATP-driven potassium channel K,ATP 

The ATP-driven potassium-selective channel (also called G-channel (8)) uses ATP to open the 
channel and to release potassium from the cell. This process is, however, also regulated by the 
presence of glucose which is a source for ATP-metabolism in mitochondria. At normal glucose 
levels the channel has a resting state open probability, which is counteracted by the activity of 
the sodium-potassium exchanger. Note that about 50% of the potassium current is mediated by 
K,ATP channels (7). At increased levels of glucose the opening of this channel is inhibited. The 
single protein current is modeled by (see also Eq.[9l): 

/k,ATP = (1 - 5'K,ATp)fi'K,ATP {V - Vk) , (14) 

with 

dgK,ATP _ g"act(7; TK.ATP^ ^K.ATp) ~ 5'K,ATP ^^^^ 

dt Tk,ATP 

Here, 7 is the glucose concentration. All other parameters are in full analogy to the voltage- 
gated sodium channel Na,V. 

Ohmic conductivity and calcium dependence: In the resting state the current is dominantly 
flowing through this channel (28) which explains why the resting potential is near of the potas- 
sium reversal potential. It is controversial whether the current-voltage relation is ohmic ( (28) 
versus (8, 9)). The latter authors show an ohmic relation in the relevant range of membrane 
potentials which justifies the present model (see (9) Fig. 2, (8) Fig. lb). The opening frequency 
is observed to depend neither on the membrane potential nor on the calcium level (9, 28) which 
justifies the present approximation Eqs.[T4]and[T5l 
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Conductivity: The conductivity of single open channels is in the order of 50pS for relevant 
membrane potentials (8) where also the dependence on the glucose level was investigated. This 
is consistent with the finding that with increasing depolarisation of the cell increasing outflow 
of potassium is observed with a single channel conductivity of (7k,atp = 54pS' (9) (this value is 
used throughout all present simulations). An approximately linear current- voltage-relationship 
is found (9). 

Activation time scale: The steady state approximation was used in two other models (4, 5). 
This points to a relatively short tk.atp- However, in an effective modelling approach (ignoring 
ATP concentrations) the time constant does not represent the time needed for the channel to 
adapt to a novel ATP-concentration but the time needed for the full metabolism process from 
glucose through ATP generation up to the changed open probability. The whole process from 
increased glucose levels up to exocytosis of insulin carrying granules take 1 minute (41), which 
defines an upper limit for tkatp- Thus, tk.atp is chosen between this time and the time scale 

Glucose-driven deactivation: As K,ATP provides the dominant part of potassium current in 
resting state it is assumed that the half-deactivation-concentration of K,ATP is around the resting 
level of glucose (7k,atp = l-2mM). The steepness of the dependence on glucose is estimated 
as kk,atp = 6mM. 

Modelling ATP-dynamics: The present simulation might be improved by inclusion of ATP- 
dynamics as done in (2, 6). The functional relationship is well supported by experiment (13). 
This will be subject of further research. 

Other ATP-driven potassium channels Other, novel, ATP-activated potassium channels have 
been observed in mice and humans, however, they differ in characteristics. They are still poorly 
defined but they are also believed to be involved in calcium oscillations. One additional channel 
found in human /3-cells has a single-channel potassium conductivity of 30^5 (28). 



5.1.6 The voltage-gated potassium channel K,V 



The voltage gated potassium channel (also called delayed rectifier potassium channel) increases 
the potassium efflux from the cell once the cell is depolarised. The potassium ions follow their 
electrochemical gradient. The model for this channel is in complete analogy to the voltage- 
gated sodium channel Eq.[9l 

Ik,y = hK,ygK,ygK^/ {V -Vk) , (16) 

with 

dt T K,V 

dhK,y _ 0'inact{V,WK,y, \K,y) — hK,v ^^^^ 
dt ^K,V 

This is largely analogous to other models of /5-cell bursting (2, 4, 5). However, the activation 
characteristics are not consistent in these models and inactivation is not considered ( (4) includes 
a different inactivating channel called /fast therein). The present model is based on single protein 
properties excusively derived from experiment. 
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Conductivity: Open channels have a single channel conductivity of gxy = lOpS (28). 

Voltage-dependent activation: K,V-channels in human /5-cells were activated with a 400ms 
depolarisations to different membrane potentials (42). Then the membrane potential was re- 
duced to -50mV and the current after this reset was monitored. The authors find a sigmoidal 
shape, centered at Vk,v = ImV and width kk,v = 8.5mV. 

Time scale of activation: The activation time constant is assumed to depend on the potential 
in (3) according to 



exp{{V - VK)/a} + exp{-(V - VK)/b} 

with c = 60ms, Vk = —75mV, a = 65mV, and b = 20mV. This function is in agreement 
with experiments where the time constant was found to vary between 8 and 37ms (42). 

Slow channel inactivation: For long depolarisations inactivation of the channels is observed. 
The characteristics are Wk,y = —25mV and width Ak,v = 7.3mV (42). However, inactivation 
is not observed within 400ms so that we must conclude that 9k,v > 400ms. The value of 
400ms is used in all simulations. 



5.1.7 Calcium-gated potassium channel K,Ca 

The calcium-gated potassium channel opens in dependence on the membrane potential and the 
intracellular calcium level. The dependence of the open probability on calcium is assumed to 
be in steady state and described by a Hill-function. Voltage-gating is described by on Ohm-like 
approach with dynamic gating (see also Eq.|9l). 

/K,Ca = 5'K,CaM7^ - "^k) ^(C,C'K,Ca,^K,Ca) , (19) 

with 

Ca Cract(^) ^4<,Ca; '^KjCa) 9K,Ca, /'OA\ 

dt rK,Ca 

All parameters are in analogy to other proteins described before. Other models that include 
detailed opening dynamics of K,Ca do not exist. The models (2,4,5) differ from the present one 
by very large Hill-coefficients (between 3 and 5) in a single Hill-function. As this is the first 
model of the K,Ca-opening dynamics it is explained in more detail. 

Ohmic conductivity: The calcium-gated potassium channel is rapidly activated upon depolar- 
isation and show a linear current-voltage relation (28). This justifies usage of an Ohm's law 
current. 

Conductivity: The conductivity of single-channels 5iK,ca = 220j9S' is rather large (28). This is 
also confirmed in rat muscle cells with values of up to 300^5* depending on the temperature (20). 

Opening dynamics: Experiments with rat muscle cells (20) show a dynamic activation in de- 
pendence on the membrane potential and of the calcium concentration (see also (28)). A sig- 
moidal opening probability is found for constant membrane potential in dependence of different 
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calcium concentrations (see (20) Fig. 6), but also for constant calcium with variable membrane 
potential (see (20) Fig. 8). 

In a first approximation the opening dynamics is modelled by a product of a sigmoidal function 
and a Hill-function as shown in Eqs.[T9] and [20l However, this approach turns out to be in 
contradiction to the data in (20). In order to reproduce these measurements more accurately, 
a dependence of the half activation calcium concentration CK,Ca(V") on the voltage has to be 
assumed (see Figure [T2|): 

dt rK,Ca 

With Cl^.AV) = expj^^l (21) 

with a = AhmV and h = 30mV, and the same time delay than for the direct voltage-dependence 
(assuming related mechanisms). This is combined with a Hill-coefficient of nK,ca = 2 and the 
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Figure 12: Voltage-dependence of half activation: The Steady-State voltage-dynamics of 
the half activation calcium concentration Cg^ca following Eq.|2l|with a = A5mV and 
b = 30mV. Note that the dynamics is relevant at physiological membrane potentials. 



sigmoidal function for the voltage-dependent activation with VK,ca = —4:0mV and ztK.ca = 
25m V. We then get the result shown in Figure [l3l A physiological interpretation of Eq.[2T]is 
not obvious. Basically this equation says that less calcium is needed to reach the same level 
of activation when the cell is more depolarised. If the dependence of the opening probability 
on voltage and calcium level are supposed to correspond to two different molecular gating- 
mechanisms, then this equation says that either depolarisation acts on both mechanisms or the 
mechanism influenced by depolarisation acts back on the calcium-dependent mechanism. A 
concise interpretation is, however, lacking, and has to be clarified in further experiments. 

Inactivation: Inactivation was only rarely observed. In rat muscle cells the channels remained 
open for several hundred milliseconds (20). Thus, inactivation of these channels is neglected in 
the model. 
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Figure 13: Opening probability of K,Ca using a product of a sigmoidal function (for 
V) and a Hill-function (for calcium). The calcium concentration of half activation is 
assumed to depend on the membrane potential according to Eq.im The results are 
shown varying calcium (left panel) or the membrane potential (right panel). There is 
good agreement with the data in Fig. 8 (20). Parameters are a = AdmV, b = 30mV, 
^K.ca = 2, VK,ca = -^OmV , and KK,ca = 25mV . 



5.1.8 Voltage-gated calcium channels Ca,L/T 

Voltage-gated calcium channels induce the inflow of calcium, following the electrochemical 
gradient of calcium, after an initial depolarisation of the cell. There are different subtypes of 
these channels which differ in their responsiveness. The dominant types in /?-cells resemble L- 
and T-type channels (28) which can be classified as HVA- and LVA-channels, respectively. The 
equations for these channels are in analogy to the voltage-gated sodium channels Eq.|9] 

/ca,L = hce.,L{l- H{C,Cc.,unce.,L))9c.,Lg^{V -Vcl) , (22) 

with 

dgCa.,L O'a.ctiV, Vca,L, «^Ca,L) — fl'Ca.L 



dt TCa.L 
dhca.,L O'inactiy, Wca,L, Aca,L) — ^Ca,L 



dt a 



(23) 



Ca,L 



For T-type channels all L's in the subscript have to be replaced by T. The term 1 — H{. . .) 
appears only for L-type channels and has to be replaced by 1 for T-type. While HVA-channels 
were included in different ways in other models, LVA-channels were neglected in all other 
models. In the following a set of experiments providing opening characteristics of L- and T- 
type channels are summarised. The simulations will rely on the single protein characteristics as 
provided in (29). 
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Single T-type (LVA) channel properties: Open T-type channels show a single protein con- 
ductivity of gca,T ~ lOpS (29) (measured in dendrites of pyramidal neurons in CAl hip- 
pocampus). They open already for membrane potentials near resting potential in the range 
of —70 to —30mV with a half opening at — 55mV (28). The threshold for opening was also 
measured in adrenal medulla endothelial cells (in whole-cell patch-mode) to be —50mV with 
Vca.,-! = — 30mV (43) (this experiment uses calcium and not barium). The threshold value of 

was also found in (29) with Vca,T = —32mV. This consistently infers Kca.,T ~ 
in order to find low activity at the activation threshold. Indeed, /tca.T = 7-OmV was found 
in (29). For inactivation Wca.,T = —Q7mV and Aca,T = Q-5mV are reported in (29). Activation 
happens on a time scale of rca,T = 10ms while inactivation is rapid (28) but is slower with 
^Ca,T = 18ms (43). 

Single L-type (HVA) channel activation: Open L-type channels show a single protein con- 
ductivity of (yfca.L ~ 27pS (29), measured in dendrites of pyramidal neurons in CAl hip- 
pocampus. In /3-cells a single-channel conductivity of 22pS was found (44). Activation is 
found between -40 and OmV ( (28), (43) in a whole-endothelial cell measurement, and (45) 
in whole-cell measurements of mouse /?-cells). The activation threshold is higher — 20my in 
CAl-pyramidal cells (29). The half maximal activation potential is Vca,L = — 18mV (43) with 
a width of Kca,L ~ lOmV. In CAl-pyramidal neurons measured with barium Vca,L = +9mV is 
found with Kca,L = 6mV (29) (however, activation depends on the barium concentration in the 
patch). Channel characteristics have also been measured in /5-cells with Vca,L = — 3.8mV^ and 
^Ca,L = 8AmV (44). The variability of whole-cell measurements of the currents is presumably 
related to different protein densities in different cell types rather that to different single-channel 
characteristics - see (43) and references in the discussion of L-type channels. Another ambi- 
guity comes from the replacement of calcium by barium for technical reasons, which induces 
differences in the measured conductivity. 

No L-type inactivation observed: Voltage-dependent inactivation is not observed for L-type 
( (28), (45) in /3-cells, (43) in endothelial cells, and (46) in bovine chromaffin cells (whole-cell 
measurement), (29) for CAl-pyramidal neurons). Thus, rca,L > Is is used. Note that the half 
inactivation potential and the steepness are irrelevant under these circumstances. Alternatively 
one could simply set /ica,L = 1- 

However, there is an experiment showing inactivation of L-type channels in cardiac cells upon 
increased calcium levels (47). Inactivation was characterised by a Hill-function with Hill- 
coefficient rica.L = 1 (thus suggesting inactivation to be mediated by the binding of calcium 
ions to a single regulating site per Ca,L-protein) and half maximal inactivation at Cca,L = 
4/iM. This behaviour is included in this model for the first time by the additional factor 
l-i7(C,Cca,L,rica,L) in Eq.[22l 

Sodium conductivity: In the same paper (44) it is also shown that the L-type calcium channels 
exhibit large conductivity of 44^5 for sodium, which is not considered in the simulations. 
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5.2 The leakage currents 

The leakage currents in Eqs.[I]and|4]are calculated in the steady state approximation. 

— -^Na = PNa,V-^Na,V + 2pNa,K-^Na,K'^Na,K + PNCX-^nCX'^NCX 
--'k = PK,ATp/k,ATP + PK,V-^K,V + PK.Ca^K.Ca - 2pNa,K/Na,K 

-</ca = PCa.L^Ca.L + PCa.x/ca.T " 2;caPNCx/NCX + PPMCa/pmcA ) (24) 

where the superscript ss denotes the steady state quantities. The calculated leakage currents 
guarantee that the cell exhibits a stable steady state and that the cell returns to its resting 
state after stimulation. In most alternative approaches the leakage currents were not consid- 
ered (e.g. (3-5)). 

The steady state currents are calculated from the steady state approximation of all kinetic equa- 
tions for the open probabilities. These are evaluated on the basis of a set of properties in the 
resting state {V^, 70, iVo, ^o, Co}- 

-^Na,K = -^Na,K (1 — H^Kq, K^a.,K, '^Na.K)) 

-^Na,V — ^act(Vo, ^Na.V, '«Na,v)5'Na,V {Vq — ^Na) 

-^NCX = -^NCX-f^(Co, C*NCX, ''^NCx) 

-^PMCA = -^PMCA-f^(C*05 C'PMCA, "^PMCa) 

-^K.ATP = (1 ~ (^actilO, 7k,ATP, I^K,ATp)) fi'K.ATP (K) — ^k) 

-^K,V = cract(V'o, 14,V, «K,v)^K?V {Vq - Vk) 

^K,Ca = ^K,Ca (^0 " Vr) H{Co, CK,Ca( Vb) , nK,Ca)cract( 14,Ca, ^K.Ca) 

-^Ca,L = <7act(K), Vca.,L, ^Ca,L)fi'Ca,L (^0 — V^Ca) 

^Ca,T = C^act(Vb, Vca,T,ft^Ca,T)^C<T (V'o - Vca) , (25) 

where CK,ca(Vo) is calculated with Eq.[2Tl The intracellular resting concentrations of ions are 
used as given in Tab.E A'o = 95mM (48), Nq = 20mM (48), Co = O.lfiM. 

5.3 Parameter lists 
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Variable 


adapted from (3) 


full model Figure[3j 


unit 


reference 


use leakage cun'ents 


no 


yes 


boolean 




use inactivation 


yes 


yes 


boolean 




PMCA feedback to V 


no 


yes 


boolean 




cell radius i?ccll 


6.1 


6.1 


fim 


(49,50) 


Cm 


10 


10 




(51,52) 


T 


310 


310 


K 




70 


1 


1 


mM 


(41) 


O'stimualtion 


none 


10 


mM 


(18,52) 


Vb 


-70 


-70 


mV 


known 




95 


95 


mM 


(48) 




5.7 


5.7 


mM 


Eq.[5] Vk = -75mV 


No 


20 


20 


mM 


(48): 20, (53): 36 




400 


400 


mM 


Eq.[5] = 80mV 


Co 


0.1 


0.1 


liM 


known 


Ce-Kt 


1.5 


1.5 


mM 


Eq.[5] W = 128mV 


use Nemst 


no 


yes 


boolean 




con'ect Vca 


no 


78 


mV 


{\l):VK = 50mV 


calcium binding sites cq 


1.0 


1.0 


mM 


(3): 0.3mM; 1.0 implies /ca ~ 0.1% 




1 


1 


fiM 


(54) 



Table 2: General framework of the simulations: The parameters as used in the full sim- 
ulation in Figure [3] and the parameters needed to reproduce the model (3) are listed. 
Note that the latter differ from the values used in (3) and are adapted to today's knowl- 
edge. References given without comment support the value used in the full model. 
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Variable 


adapted from (3) 


full model Figurel3l 


unit 


reference 




none 


0.00003 




(24) 




none 


33.3 


mM 


(22) 


77, TS.T ^ Xi' 

' '^iN a,l\. 


none 


2 


# 


(22) 


-^Na,K 




20 




C221 






2 


# 


(22) 


^Na.K 




1.5 


Na/K 


known. 


5k,atp 




54 


po 


(8 9 28) 


'TCjATP 




I 




vfl ntiA'fl^/ clou/ 


TK ATP 


none 


1.2 


mM 


lai'ge cuiTent at resting state (7) 


kk,atp 


none 


6 


mM 


estimated 


9K,V 


10 


10 


pS 


(28) 


tk.v 


Eq.[T8] 


Eq.[T8] 


s 


(42) 


half max tk v 


0.03 


0.03 


s 


8ms < TK,v < 37ms (42) 


Vk.v 


-15 


1.0 


mV 


(42) 


kk,v 


5.6 


8.5 


mV 


(42) 




none 


400 


ms 


(42) < 400ms 


VFk,v 


none 


-25 


mV 


(42) 


Ak,v 


none 


7.3 


mV 


(42) 


9k, Ca 


220 


220 


pS 


(28), (20): /(T) 


CK,Ca 


0.001 


Eq.|2l] 


mM 


(20): /(V) 


"K,Ca 


3 


2 


# 


(20) 


TK.Ca 


none 


100 


ms 


estimated 


VK.Ca 


none 


-40 


mV 


(20) 


KK,Ca 


none 


25 


mV 


(20) 


9Na,V 


none 


14 


pS 


(25) 


ma,V 


none 


Eq.ni] 


s 


(11,19) 


14*a,V 


none 


-35 


mV 


(26) 


tNa,V 


none 


8.0 


mV 


(26) 


^'Na.V 


none 


4.6 


ms 


(11) 


WNa.V 


none 


-100 


mV 


(19,26,28) 


ANa,V 


none 


20 


mV 


(19,26,27) 



Table 3: Properties of single sodium and potassium conducting transmembrane proteins: 
Properties of sodium and potassium conducting membrane proteins. References given 
without comment support the value used in the full model. A comment of the form /(X) 
denotes that the reference claim a dependence on the variable X. 
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Variable 


adapted from (3) 


full model Figure [3j 


unit 


reference 


^NCX 


none 


-0.5 


fA 


(30) 


Cncx 


none 


1.8 




(31) 


riNCX 


none 


1 


# 


(10) 




none 


3 


Na/Ca 


(31) 


^PMCA 


0.01 


0.01 




(30) 


CpMCA 


1.8 


0.1 




(40) 


^PMCA 


2 


2 


# 


(38) 


'^PMCA 


1 


1 


Ca/ATP 


(37) 




27 


27 


vS 


(29, 44) 




1 


6 


Tfl.S 


(29): fast 


Vca L 


0.0 


0.0 




(29,43,44) 




12 


12 


rnV 


(28,29,44) 


6Ca,L 


10 


10 


s 


(29): no inact. 


W'Ca,L 


100 


100 


mV 


(29): no inact. 


-\Ca,L 


10 


10 


mV 


(29): no inact. 


Cca.L 


none 


4 




(47) 


'^Ca,L 


none 


1 


# 


(47) 


9Ca,T 


none 


10 


pS 


(29) 


TCa.T 


none 


10 


ms 


(43) 


Vca,T 


none 


-30 


mV 


(29) 


«Ca,T 


none 


7.0 


mV 


(29) 


^Ca,T 


none 


18 


ms 


(43) 


VKca.T 


none 


-67 


mV 


(29) 


-^Ca,T 


none 


6.5 


mV 


(29) 



Table 4: Properties of single calcium conducting transmembrane proteins: Properties 
of calcium conducting transmembrane proteins. References given without comment 
support the value used in the full model. 



Variable 


adapted from (3) 


full model Figure [3j 


unit 


Pk,ATP 


none 


0.092 


/ ^m'^ 


Pk,v 


0.24 


8.0 


/ ^m? 


PK,Ca 


0.08 


0.45 


/ ^m? 


PNa.K 


none 


2000 


/ ^m? 


PNa.V 


none 


1.15 


/ ^m? 


Pncx 


none 


7.5 


/ ^m'^ 


Ppmca 


10000 


1350 


/ lim? 


PCa.L 


0.1 


0.9 




PCa,T 


none 


0.1 





Table 5: Densities of ion-conducting transmembrane proteins: The densities of trans- 
membrane proteins as used in the two simulation set-ups are listed. These parame- 
ters were used as fit parameters for the simulations. Thus, no references are given. 
Note that the large PIVICA-density is a result of an unrealistic half activation calcium 
concentration used in (3) (see also Tab.H]). 
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